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ABSTRACT 


Theoretical investigations concerned with efficiency of rf power transfer to 
a cylindrical plasma are presented. The relationships between the six magnetic and 
electric wave field components are derived for the quasistatic case. Through suitable 
approximations, relations between the various wave field components are developed 
so that each of the components may be expressed in terms of quantities measurable 
in the laboratory. A method is proposed for determining wave power flow from 
magnetic probe measurements. 

Experimental measurements of rf power coupling to the plasma are described. 

A measurement of phase shift is described, with initial results indicating very large 
\ II o Results from the wavelength measurement are compared to theoretical calcu- 
lations of \if based on measured electron densities. Observations of complex fine 
structure in the wave magnetic signal are also reported. 

Evidence of non-resonant coupling of rf power to the plasma is discussed. 

The so-called "Mode X" which is thought to be responsible is described, along with 
suggestions for its elimination, involving improved Faraday shielding of the Stix Coil. 

An investigation of waves which propagate at low magnetic fields (B 0 1000 G) 

is reported. These waves, which have been previously referred to as the so-called 
"fast" or compressional hydromagnetic waves, are found to propagate in "bands", 
or ranges of magnetic field. These "propagation bands" are separated by "stop bands", 
or ranges of magnetic field for which the waves appear attenuated. 

It is found that during the strongest wave propagation, in the central 
''propagation band", electron density increases by more than a factor of two. Plasma 
diamagnetism is found to decrease, however, indicating a net loss of energy per ion. 

An optimum pressure is found, at which wave amplitude and electron density are at a 
maximum. 

The study of nonlinear frequency phenomena is continued and extended to 
low magnetic field waves. Distortion of magnetic signal is found to increase strongly 
during wave attenuation. A comparison of the frequency spectra of wave magnetic 
signal ib) and Stix Coil current during the rf heating cycle verify that the harmonic 
frequency generation occurs in the plasma rather than in the rf transmitter. 
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PROPAGATION AND DISPERSION OF HYDROMAGNETIC AND ION CYCLOTRON 
WAVES IN PLASMAS IMMERSED IN MAGNETIC FIELDS 

Jimmy G. Melton and Nathan B. Dodge (Prof. Arwin A. Dougal) 

I . INTRODUCTION 

Damping of the ion cyclotron wave in a magnetic beach offers a promising 
method for heating plasma ions to extremely high temperatures. Since the process 
heats predominantly the ion constituents of the plasma, it can be used to heat the 
plasma ions directly, without relying on energy transfer through electron-ion 
collisions to achieve high ion temperatures. 

1 2 

It has been predicted and verified that power from an rf oscillator can 
be coupled to the ion cyclotron wave in a plasma with over 50% efficiency, the 
losses being primarily copper losses in the coupling arrangement. It has further 
been observed that about 90% of the energy in the wave can be transferred to the 

3 

ions in the damping region. Thus, the overall transfer of energy from the driving 
circuitry to ion energy can be made as high as 45% effective. 

The primary objective of ion cyclotron resonance heating studies is to be 
able to heat the overall plasma as hot as possible with the power available. This 
involves optimizing each stage of the heating process with regard to the plasma 
temperature attained. The heating process may be regarded as having three stages: 
(1) coupling of energy to the plasma column, (2) transfer of the energy from the 
coupling region to the damping region via the ion cyclotron wave, and (3) the 
damping of the wave and the resulting ion heating. 

A complete study of optimizing the heating process must include all of 
the losses and limitations involved. This includes energy lost due to the loss 
of particles, by anomalous diffusion or otherwise, energy lost due to undamped 
waves, or lost due to coupling to undesired waves, or heating limitations due to 
finite temperature effects. However, much is still to be learned about the energy 
which eventually heats the plasma ions . 


1 
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The studies in progress, and described in this report, are concerned with 
determining heating efficiency of the ion cyclotron wave. Section II of the report 
presents analytical studies of the electromagnetic wave fields in the quasistatic 
approximation, with the objective of establishing means to measure experimentally 
the wave energy in the plasma wave. Also in Sec. II are experimental measure- 
ments performed to determine the rf power being coupled into the plasma. It has 
been found that little of the available power is being coupled into the plasma. 

Evidence is presented which relates the inefficient coupling to the "Mode X" 

4 

coupling phenomena observed and effectively eliminated at other facilities. 

Section III presents the results from an experimental investigation of the 
waves which propagate below the ion cyclotron magnetic field. One of these 
waves has tentatively been identified as the "fast" hydromagnetic wave. It has 
been found that rf power can be coupled to this wave with about 10 times the 

5 

efficiency as to the ion cyclotron wave. It has also been found that the fast 
wave can be used to heat the plasma. This offers possibilities of preheating the 
plasma with the fast wave or of heating with the fast wave and the ion cyclotron 
waves simultaneously. 

It was found that low magnetic field wave propagation occurs in "bands," 
or for separated ranges of magnetic field. It was also found that wave propagation 
and increased density occur together. This suggests that propagation of the wave 
at low fields increases the ionization. 

Also presented in Sec. Ill are further observations on the generation of 
frequency harmonics when waves damp in the magnetic beach. It was reported 

g 

in the previous semiannual report that generation of harmonics of the wave 
frequency was found to be associated with wave propagation and damping in the 
magnetic beach region. Further investigations have revealed that a significant 
fraction of the wave energy may be converted to new modes which escape cyclo- 
tron damping. 
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Investigations are being continued into the flow of rf power in ion cyclo- 
tron wave propagation and damping. It is through studies such as these that 
the evaluation of the ultimate capability for plasma heating by ion cyclotron 
resonance absorption may be accomplished. 
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DETERMINATION OF POWER FLOW IN ION CYCLOTRON RESONANCE 
PLASMA HEATING, Jimmy G. Melton and Nathan B. Dodge (Professor 
Arwin A. Dougal) 
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The central objective in plasma heating by ion cyclotron resonance heating 
is to maximize the amount of power being used to heat the plasma. In present- 
day experiments, the objective has become rather to optimize the plasma heating 
power. The problem becomes in practice twofold: to optimize the plasma heating 
process and to optimize the delivery of power to the heating process. Having 
achieved these, it is hoped that the overall process will be efficient enough to 
enable heating a large volume of plasma to thermonuclear temperature. 

The ion cyclotron resonance heating process proceeds as follows: energy 
is coupled to the plasma column from a radio frequency source operating at or 
near the ion cyclotron frequency of the plasma. Some of this energy goes into 
the ion cyclotron wave, which propagates for 6D< Sh' , where Ac is the 
ion cyclotron frequency, Al ~ * T ^ e ion cyclotron wave is then made 
to damp via the ion cyclotron damping process by making d) = At . As the 
wave damps, it delivers its energy to the ions present, which then thermalize 
through collisions . 

From this brief description it can be seen that there are three phases to 
be considered in optimizing plasma heating; (1) the coupling of energy into the 
ion cyclotron wave, (2) the delivery of energy by the wave to the damping region, 
and (3) the transfer of energy from the wave to the ions of the plasma. Each of 
these three phases must be adjusted to yield maximum plasma heating for the 
input power available. 

Figure II- 1 is a simplified block diagram of power flow. It shows only 
the dominant features of power flow, omitting many of the details such as elec- 
tron heating. The power from the rf source goes into ionizing and heating the 
plasma, into ohmic losses in the coupling device, and into the plasma. Of that 
part which goes into the plasma, the percentage going into ion cyclotron waves 
can be made quite large by suitably adjusting the experimental conditions. The 
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power going into the plasma can be determined experimentally from measurements 
of incident and reflected power and from tuning conditions . In some cases the 
percentage of power going into ion cyclotron waves can also be determined by 
measurements of the magnetic signal associated with the wave. Similarly, 
measurements of the plasma diamagnetism can be used to determine plasma 
heating and temperature. 

From measurements of this type , it is hoped that the efficiency of power 
transfer of each phase of the process can be determined. By comparing such 
measurements with results obtained from theoretical analysis, the optimum condi- 
tions for plasma heating can be determined, which would then afford an overall 
evaluation of the ion cyclotron resonance heating method. 

Such optimization studies can be carried out using medium-powered 
experiments with experimental arrangements which stress flexibility. The results 
thus obtained can then be tried on high-power experiments of less flexible 
design. 

It is to this objective of optimizing each phase of power transfer that 
these studies are directed. Experimental techniques are being developed to 
measure the power at each step in the energy transfer process. Analytical studies 
are being conducted in support of these techniques. 

Section II-A presents a derivation of the relations between the electric 
and magnetic field components for ion cyclotron waves in the quasistatic approxi- 
mation. With these relations, one can theoretically determine the wave energy 
density and power flow of ion cyclotron waves from measurements on the elec- 
tric and magnetic fields of such waves . 

A. Quasistatic Wave Field Amplitudes 

At the present time, the wave coupling phase of ion cyclotron resonance 

1 2 

heating experiments is relatively well understood. Several analyses ' of the 
wave coupling process, in different approximations and varying degree of 
detail, have found fairly good agreement with experiment. There remains. 
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however, one aspect of experimental measurements on coupling efficiency which 
needs to be performed with more accuracy. Since not all of the power going into 
the plasma goes into the ion cyclotron wave, it is necessary to measure the 
energy flux of the wave in order to know how much was coupled into the plasma 
in some form other than ion cyclotron waves. Such measurements of the energy 
in the wave require quantitative measurements of the wave field amplitudes and 
a knowledge of the local plasma conditions for the measurement position. This 
section presents the relationships which exist between all six of the wave mag- 
netic and electric field components in the quasistatic approximation (which seems 
to be valid for the conditions of the laboratory plasma under consideration). 


The electric field of a wave in the plasma is given by the solutions of 
the wave equation 

V X (vxE) + ~r K'E — ° 


which satisfy the boundary conditions. A similar relation holds for the magnetic 
wave field B . Here K is the effective dielectric tensor of the plasma. 


For a cylindrical cold plasma of radius a and dielectric constant 

-n 3 

the field parallel to a constant confining field B q £, are given by. 
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The boundary condition equations are solved simultaneously with the dispersion 
relation 


7' 




(II- 4) 


V 2 ) 2 

to obtain v > and ^ . However, to solve the problem in its general 

form requires numerical computation. 


The problem can be simplified considerably by suitable approximation to 

yield a dispersion relation with only one J solution. In the quasistatic 
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approximation, in which we are interested, 4/ is large, and 
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In this approximation, may be neglected so that, by taking — O , 

which is important only near the boundaries, one can derive the following field 
relations , 
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The above relations are sufficient to determine each of the wave field com- 
ponents in terms of theoretically measurable quantities and other field components. 
However, it is also possible to express each field component in terms of one of 
the others; for instance. 




E, 




Ov - a 
^ W 


(II- 6) 


Expressed in this way, within the limits of the approximation, i.e. ^ 1 / 

it is seen that E is the dominant electric component, with both E and E 
r z y 

being much less than E r , whereas Bg>^ > B f holds for the magnetic field 
components . 


Having expressed the wave field components in this form, it is now 
possible to propose a method of determining wave power flow from magnetic probe 
measurements of the wave. The method depends on the quasistatic approximation 
being valid (though if it were not valid, wave field relations could still be 
derived which would be different from those presented above) , that the pertinent 
plasma parameters are known (density and magnetic field), and, most importantly, 
that no more than one transverse mode exists and the wavelength parallel to 
the magnetic field can be determined. 

If these conditions are met, then the transverse wavenumbers ^ may 
be determined from the relation, 

B Ve* = " (II - 7) 

Knowing /0 , the wave field relations can be used in Poynting's theorem for 

the plasma, as presented in Ref. 4, to determine the wave energy content for 
quasistatic ion cyclotron waves . If the quasistatic condition is not a good 
approximation for the plasma, a similar analysis will yield the appropriate field 
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relation to be used. The critical assumption, however, is that the wave is com- 
posed of a single transverse mode. If this is not the case, other means must be 
used to determine wave energy density. 

B. Experimental Measurements of Power Coupling 

The experiment performed on power coupling to the plasma consisted of 
measuring incident and reflected power to the plasma, together with simultaneous 
measurements of the current flowing in the Stix coil, first with no gas and then 
with plasma in the system. Analysis of the data reveals the power being lost to 
the resistance of the coil, and the combined power going into the coil resistance 
and into the plasma. 

The experimental arrangement is shown in Fig. II— 2 . It consists of a 
Pyrex discharge tube 6.25 cm in outside diameter located in a solenoidal magnet 
system. The magnetic field is adjustable up to 10 kG. The field configuration 
has modest 10% mirrors on each end, a 40 cm region where the field is uniform 
to +1/2%, and a "magnetic beach" of 50 cm length, which has a minimum field 
of 82% of that in the uniform region. The distance between mirror peaks is 130 
cm. The magnetic field configuration is shown in Fig. II— 3 . 

Preionization of the deuterium gas is provided by either or both of two rf 
systems; (l)a40kW oscillator at 7.5 MHz capacitively coupled to the plasma 
through two copper straps at the mirror peaks , and (2) a 1.5 kW rf amplifier at 
21.3 MHz connected between internal electrodes located at the ends of the system 
outside the mirrors. 

The main rf unit is a 30 kW pulsed amplifier at 5 . 8 MHz which is con- 
nected through directional couplers and a tuning box to a Stix coil surrounding 
the discharge tube in the uniform region of the magnetic field. A Rogowsky belt 
for measuring coil current is located around the feed strap to the Stix coil. 

The Stix coil is center-fed, with four phase-reversed sections of three 
turns each. It has a wavelength of 16 cm. Electrostatic shielding is provided 
by an ungrounded "squirrel cage" Faraday shield, consisting of a number of axial 
strips of 0.060" copper 3/8" wide bonded between mylar sheets.'* 



WATER COOLED MAGNET COILS 



Fig. II— 2 Block Diagram of ICRH Experiment 
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On the beach end of the discharge tube, a movable probe carriage 

provides access to the plasma for magnetic probes by means of a re-entry glass 

tube--bellows arrangement. Also located in the beach region is a diamagnetic 

loop around the outside of the discharge tube to provide a measure of the 

plasma density- transverse temperature product. At the time this experiment was 

performed the diamagnetic probe was not calibrated and was used only to provide 

a relative measure of plasma heating. Similarly, an 8-mm microwave inter- 

11 -3 11-3 

ferometer has been used to measure densities from 2x10 cm to 8 x 10 cm , 
although density measurements were not performed at the same time the power 
input measurements were made. 

The input power experiment consisted of simultaneously recording incident 

power, reflected power, and Stix coil current, first with no gas in the system, 

and then with plasma. Magnetic probe signal (b ) and diamagentic probe signal 

z 

were monitored to verify wave propagation and plasma heating. 

The effective plasma resistance can be calculated from power and current 
measurements for the vacuum and plasma cases by means of the relation 

R = R c + R f = (P;-P^)/le (n-8) 

where R is the Stix coil resistance, R is the effective plasma resistance, P. is 
c p 1 

incident power, P is reflected power, and I is the Stix coil current. (For 
r c 

further details see Ref. 5.) 

Figure II— 4 shows the data obtained. For these oscillographs, the pre- 
ionization was not used, such that the plasma was created in the middle of the 
main ICRH pulse. Data obtained are plotted in Fig. II— 5 . 

The data reveal that only a small percentage (less than 10%) of the power 
is going into the plasma. This point is considered later in the report, with a 
possible explanation and remedy. 
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A. Wave Magnetic Signal 
0.2 v/cm 

B. Stix Coil Current 
about 100 amp/cm 

C. Reflected Power 
1 v/cm 

D. Incident Power 
1 v/cm 

Gas Pressure: 0.4 millitorr hydrogen 
Magnetic Field: 6250 gauss 




Fig. II— 4 . Oscillographs of Incident and Reflected Power 
and Stix Coil Current. 

The time scale on these oscillographs is 200 p s/cm. 
Before 400 psec, there is little or no plasma present. 
Breakdown occurs at 400 psec. 
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FIG. E-5 POWER COUPLING TO PLASMA vs MAGNETIC FIELD 
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In order to determine how much of the input power was appearing in the ion 

cyclotron wave, a magnetic probe (b ) was used to measure wave magnetic signal 

z 

amplitude and to determine the axial wavelength for the waves . The wavelength 
measurements were performed by measuring phase shift of the wave as the 
magnetic probe was moved through the beach region. To measure phase shift, 
the Stix coil current, measured by the Rogowsky belt, was used as a phase 
reference, and displayed on the same dual beam oscilloscope as the wave magne- 
tic signal. Typical data are shown in Fig. II— 6 . 


No phase shift was observed for the range of conditions investigated. It 
is estimated that the axial wavelength would have to be greater than 200 cm 
before no phase shift could be detected by this method. Aj| = 200 cm implies a 
of 7T x 10 2 or ^ZZIO ^cm 1 . To determine if the phase measurements 
are possible, the curve of Fig. II-7 was plotted, using the cold plasma dispersion 
relation for the ion cyclotron wave. 
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A value of j^> =1.0 corresponds roughly to the lowest order radial mode. For 


Ts?j_ - 1.0, and a density of greater than lO^cm ^ it is seen that ~^ , li 


is on the 


order of 10 or the axial wavelength is 20 cm, 


In order for the wavelength to be as long as the phase shift measurements 

11 -3 

indicate, either the density must be much less than 10 cm , or the magnetic 
probe is responding either to noise or to a long wavelength "standing -wave" 
oscillation. Since it is not yet settled which situation exists, further investiga- 
tions will be conducted to resolve the question of wavelength. 


During the experiment to determine wavelength, it was observed that the 

wave magnetic signal (b ) can possess a very complex fine structure. This fine 

z 

structure at times became quite pronounced. Typical fine structures are shown 
in Fig. II- 8 on three different time scales. There are oscillations at 25 MHz, 
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Phase- Shift Measurements 

Upper trace: Wave magnetic signal 
0 . 5 v/ cm . 

Lower trace: Stix Coil Current 
@ 50 amp/cm. 

Time scale: O.ljasec/cm. 

Conditions: p = 0.3 mT of D 
Bq = 4100 gauss 

Top pair: Magnetic probe located 
at Z = 51 cm . 

Middle pair: Z = 41 cm. 

Bottom pair: Z = 31 cm. 


nd 

Fig. II— 6 * Determination of Axial Wavelength of 2 
Harmonic Ion Cyclotron Waves. 

Wavelength is measured by determining the 
phase-shift between the current in the Stix 
Coil and the magnetic probe b^. 

The probe position is indicated by its distance 
Z from the center of the Stix Coil. 

The axial wavelength of the Stix Coil is 16 cm. 
No phase shift through the plasma is evident 
in these oscillographs. 



kj. (cm~2) 


Fig. II- 


7 . Axial Wave Number vs Transverse 
Wave Number for Cold Plasma Ion 
Cyclotron Wave , with Electron 
Inertia Effects . 
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Wave Magnetic Probe Signal 
(1 .) 0.2 psec/cm 

(2.) 2.0jisec/cm 

(3.) 20 /asec/cm 

Conditions: p = 0.3 mT in 

B = 4500 gauss 
Z°= 45 cm 


Fig. II- 8. Fine Structure on Wave Magnetic Signal. 

Showing the shape of the wave magnetic signal 
on successively slower time sweeps for conditions 
which give a propagating wave. 

Probe is positioned 45 cm from the center of the 
Stix Coil, in the magnetic beach region. 

The top oscillograph shows oscillations at the 
driving frequency, 5.8 MHz, and at about 25 MHz. 

The middle oscillograph shows oscillations at 
about 280 KHz . 

The bottom oscillograph shows oscillations at 
about 30 KHz. 
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5.8 MHz (the driving frequency), 27 8 KHz, and about 30 KHz. Similar oscilla- 

6 7 

tions have been observed in other experiments and previously on this experiment. 
The cause and the nature of these oscillations is not known at present. It is 
interesting to note that the oscillations seem to be of the relaxation type. 

C . Evidence of Non-resonant Coupling of Power to the Plasma 

In the previous section it was noted that the power coupled to the plasma 

was lower than expected. Of greater significance is the observation that the 

power to the plasma does not peak near B = B . This suggests that the power is 

c 

coupling to the plasma by some mechanism other than to the ion cyclotron wave. 
Reports on experiments at other facilities have mentioned similar non- 

g 

resonant coupling phenomena. On the Model C Stellarator a similar occurrence 
of non-resonant coupling was remedied by improving the shielding of the E 

z 

component Stix coil. The mechanism which has been advanced for "Mode X," 
as it is called, is that the fringing E fields into the plasma accelerate the 

Z 

electrons along the magnetic field lines, so that the resulting space charge field 
drives the ions into the wall of the discharge tube. The significant improvement 
obtained by introducing a Faraday shield between the plasma and the Stix coil 
lends support to this mechanism. 

In addition to the power input measurements, two observations suggest 
that "Mode X" is limiting the power to plasma waves in the present studies. 

The first is that plasma density measured by an 8-mm microwave interferometer 
is present even when no gas is in the system. The second is that the current 
drawn by the ionization vacuum pump increases when the system is pulsed several 
times in rapid succession. Both of these observations suggest there is a high 
particle influx rate from the tube wall during the ICRH pulse. 

The apparent solution is to improving the present shielding arrangement 
of the Stix coil. An improved Faraday- shielded Stix coil has been designed and 
constructed and will be installed on the experiment. With the new coil in place, 
studies will be continued into the efficiency of coupling power into the ion 
cyclotron wave . 



22 


References for Section II 

1. T. H. Stix, "Generation and Thermalization of Plasma Waves , " Phys . 

Fluids , l, 4, p. 308 (1958). 

2. R. R. Woollett, "Energy Addition to an Atomic Hydrogen Plasma at Off- 
Resonant Conditions," NASA Technical Note TN D-2982, September, 1965. 

3. D. G. Swanson, R. W. Gould, and R. H. Hertel, "Experimental Study of 
Compressional Hydromagnetic Waves , " Phys. Fluids ]_, 2, p. 269 (1964). 

4. J. G. Melton and A. A. Dougal, "Theoretical Investigations of Quasi- 
static Ion Waves and Wave Energy Densities," S.S.R. No. 8 on NASA Grant 
NsG-353 , pp. 34-50, January 15, 1967. 

5. M. Kristiansen, J. G. Melton, F. C. Harris, N. B. Dodge and A. A. 
Dougal, "Experimental Investigation of Wave Coupling Structure," S.S.R. 
No. 6 on NASA Grant NsG-353, pp. 36-51 and pp. 20-26, January 15, 

1966. 

6. R. S. Pease, S. Yoshikawa and H. P. Eubank, "The Confinement of Plasma 
Heated by Ion Cyclotron Resonance in the C Stellarator , " MATT-Q-23, 

p. 7, April, 1966. 

7. M. Kristiansen, F. C. Harris, N. B. Dodge, J. G. Melton and A. A. 

Dougal , S.S.R. No . 5 on NASA Grant NsG-353, pp. 35-44, July 15 , 1965 . 

8. S. Yoshikawa, M. A. Rothman and R. M. Sinclair, "Nonresonant Coupling 
of ICRH Power to C-Stellarator Plasmas , " MATT-Q-22, pp. 28-35, April, 
1965. 


23 


III, EXPERIMENTAL INVESTIGATION OF HYDROMAGNETIC WAVE PROPAGATION 
AT LOW MAGNETIC FIELDS, Nathan B. Dodge and Jimmy G. Melton 
(Professor Arwin A. Dougal) 

Observations of wave propagation at very low magnetic fields have pre- 
viously been reported'*' from this facility. The waves which have been seen at 
low magnetic fields (B ^ 1000G) have been referred to as the so-called "fast" 

hydromagnetic wave.*' Since heating and enhanced ionization of plasma may 

2 3 

occur during the "fast" wave propagation, ' and since all hydromagnetic wave 
phenomena are of interest in the study of radio- frequency plasma heating, a study 
has been made of these low magnetic field wave phenomena, the initial results 
of which are reported here . 

Part A of this section describes the behavior of the wave phenomena, 
including onset and attenuation characteristics , and sensitivity of the wave to 
neutral pressure and magnetic field. Plasma heating and electron density during 
the wave propagation are discussed in Part B. Part C extends the nonlinear 
frequency phenomena investigation initiated in the previous report period to the 
study of harmonic generation in low magnetic field waves . Part D consists of a 
discussion and interpretation of the three sections preceding, and Part E includes 
plans for completing the low field investigation and continuing the study of non- 
linear plasma phenomena. 

A. Hydromagnetic Wave Propagation at Low Magnetic Fields 

At lower magnetic fields such that o»IL (B£l.O kG) wave phenomena 

1 c 

have been observed and reported. These phenomena were taken to be due to 
the so-called "fast" hydromagnetic wave . Recent studies, considered here, 
confirm the existence of extensive wave propagation in one or more bands. Data 
presented here is confined to hydrogen. 

The wave propagation depends rather sensitively on neutral gas pressure. 

At 0.3 mTorr, as seen in Fig. Ill- 1 , there is a gradual decay in wave magnetic 

signal, b , as the static magnetic confining field B is reduced. Low- field 
z o 

propagation (B 1000 G) has not been observed below 0.3 mTorr. 



Fig. Ill- 1 . Wave Magnetic Signal (b 
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Onset of propagation appears to occur normally at about 0.4 mTorr, around 

500 G, as shown in Fig. III-2a. For comparison, b at 0.3 mTorr, 500 G is also 

. z 

shown in Fig. Ill— 2 b . Wave magnetic signal (b ) is ragged in the former figure, 

z 

indicating that the conditions for propagation are marginal, at best. The wave 
appears to attenuate at 400 G, but again "tries" to propagate at 300 G, as in 
Fig. Ill- 2c. There appears to be marginal propagation down to around 200G, but 
below this field the wave field gradually decreases until the wave appears to be 
cut off around 75 G, (Fig. Ill- 2d). 

The appearance of two fuzzily separated propagation bands is further 

illustrated by data taken at 0.5 mTorr. A plot of magnetic signal versus confining 

field reveals two fairly broad regions of propagation separated by a narrow, 

indistinct "stop-band" or region of signal attenuation, as shown in Fig. Ill— 3 . 

Again, there is a large amount of variation in the amplitude of b , as indicated on 

z 

the graph, in the propagation bands, so that evidently conditions are still near- 
marginal for propagation. 

It is interesting that the "stop-band" appears to have narrowed somewhat 
from the band at 0.4 mTorr and is only about 50 to 100 G wide. The stop bands 
and propagation band retain the same configuration at 0.6 mTorr, but the stop 
band appears to be wider in this case, and the propagation band around 500 G 
seems somewhat wider. The trend continues at 0.7 mTorr. 

Figure III— 4 shows b^ versus B q at 0.8 mTorr. The magnetic signal 
shows a broad peak around 500 G and two lesser peaks at 1 kG and 75 G, with 
broad stop bands of around 200 G between each one. There now appear to be 
three regions of propagation: (1) a weaker band below 100 G, (2) a broad, rather 
sharply defined band from 300 to about 650 G characterized by high amplitudes 
in b , and (3) a band beginning at around 800-900 G. (Although this study was 
not carried out at magnetic fields above 1.0 kG, previous work indicates this 
band extends upward to about 2.2 kG at 0.8-0. 9 mTorr in hydrogen gas.) 

As implied by the above description, the stop bands extend from 100 to 
300 G and from about 650 to 900 G. Studies so far have extended to 1 . 2 mTorr in 




a. (iower trace) Onset of propagation, first band. 

p = 0 . 5 mT Z = 50 cm 

gas Oscilloscope Sweep: 200 yusec/cm 

B = 500 G 
o 



b. b^ (lower trace) No propagation, 
p = 0 . 3 mT Z = 50 cm 

H 2 gas Oscilloscope Sweep: 200 yusec/cm 

B = 500 G 
o 
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c. b (lower trace) Onset of propagation, second band, 
z 

p = 0 . 4 mT Z = 50 cm 

H 2 gas Oscilloscope Sweep: 200 yusec/cm 

B = 300 G 
o 



d. b^ (lower trace) Strong attenuation, 
p = 0 . 4 mT Z = 50 cm 

H 2 gas Oscilloscope Sweep: 200 yusec/cm 

B = 75 G 
o 


III- 2 . b^ During Propagation and Damping. 

(Z = Distance to Center of Stix Coil) 
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hydrogen. Above 0.8 mTorr, there is a gradual disappearance of the lowest 
propagation band, while the central band broadens to include a range of 300- 
700 G at 1 . 2 mTorr. The lowest band is no longer in evidence at this highest 
pressure so far studied. The higher band is still in evidence, although it is 
about 500 G wide, extending from 800 or 900 G to 1.4 kG, indicating that a 
considerable narrowing of the band has taken place. (The upper limit of the 
band was verified from a previous study, as mentioned above for (3).) 

Further data will be taken above 1.2 mTorr. The studies referred to above 

indicate, however, that above 1 or 2 mTorr a "flattening" of the peaks in the 

b^ - B q curves and a "filling" of the valleys occurs. By about 10 mTorr, the b^ 

signal is more or less level over the range of B between 50 G and 1000 G. It 

o 

is assumed that the reason for the gradual "smearing out" of wave phenomena 
above 1-2 mTorr is due to ion-neutral collisional effects, about which more will 
be said in Part D. At present, no investigations can be made below 50 G as this 
is the lower limit of the dc power supply. 

B. Electron Density and Diamagnetic Signal Measurements of Low- Field 
Wave Phenomena 

4 

As has been reported previously, a microwave interferometer has been 
installed which allows precise measurements of electron density. It is therefore 
possible to accurately monitor the percentage ionization in the laboratory plasma 

during the radio frequency heating cycle. With the provisions for measurement 

5 6 

of diamagnetic signal already available, ' it is a simple and straightforward 
procedure to calculate the relative change in ion temperature. Recall^ 


Te x +\ =(&) 


V, RC 
NA % m 0 k - 


(HI-1) 


where and~Tij_ are the perpendicular electron and ion temperatures, respec- 

tively, n is the electron or ion density (assuming that n = n. = n), V is the 

e i s 

integrated diamagnetic signal, the RC product represents the time constant of the 
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integrator, B q is the static magnetic confining field, N is the number of turns 

in the diamagnetic probe coil, A is the cross-sectional area of the plasma 

column, £ is the shape factor of the cross section, ^Cio is the magnetic per- 

_ 7 

meability ( = 4*Tr x 10 hy/m ), and X is Boltzmann's constant ( = 1,38 x 
10~ 23 J/°K). 


1 


It may be seen that the integrated diamagnetic signal V g is proportional 


to Hr - ) by rewriting Eq . (Ill- 1 ) , 

D 

o 


Vo = 


HA£ 

■RC 




(i) 


(HI- 2) 


Thus, as B q is lowered, as in this experiment, V g would normally increase, 

even though the n(T + T. ) product remained constant. In order that the 

e X U- 

variation of diamagnetic signal be strictly due to a change in n (ionization) or 
T and T. (particle energy), the amplitude of V was corrected by multiplying 

G pL \ - S 

each value of V by the corresponding value of B in kG, so that the term "diamag- 
s o 

netic signal" will subsequently refer to the V B product, which varies only with 

s o 

density and temperature. 


Since the microwave interferometer operates at very low power levels, and 
since the microwave equipment is operated in close proximity to the plasma 
environment, a large amount of radio noise contamination of the microwave signals 
was experienced. In order that difficulties in interpreting the composite signals 
would not arise the klystron reflector voltage was square-wave modulated so that 
the klystron in the microwave system was pulsed on and off, as seen in Fig. Ill— 5 . 
The difference in the "step" was measured to register the phase shift and therefore 
the plasma density . 

Electron density is plotted versus B q at 0.3 mTorr in Fig. Ill- 6. It is 

interesting that the density falls very nearly as b , and it was found that the 

• z 

correspondence between b and electron density was close in all the cases of 

z • 

interest. Figure III— 7 depicts b and electron density plotted versus B at 

z o 



p - 1.0 mT Z = 40 cm 

H 2 gas Oscilloscope Sweep: 200 yusec/cm 

B = 300 G 
o 


Fig. Ill- 5. Typical Display of Microwave 

Interferometer with Square Wave 
Modulation of Reflector Voltage 
to Compensate for Noise Level. 




Fig. Ill- 6 . Electron Density at 0.3 mT. 
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0,8 mTorr (as mentioned earlier, all data mentioned in Section III was obtained 
for hydrogen gas) . 

♦ 

Note that there is a close correspondence between the peaks in b and 

z 

those in electron density. There is an especially close correpondence between 

the peak in b at around 500G (which is the center of the middle "propagation 
z 

band") and the peak in electron density. This was also true at 0.5 mTorr and 

1.0 mTorr, although the correspondence is not as close in the lower propagation 

band (50-100G). During strong propagation at 0.8 mTorr when the magnetic 

signal is almost triple the "stop-band" values (at B & 500G), the electron 

11 -3° 12 -3 

density has increased from a low of 5.5 x 10 cm to 1.4 x 10 cm , which 
corresponds to an increase from about 1 % to more than 2 % in total ionization. 
Similar increases are evident at other neutral gas pressures for the propagation 
band centered at 500G. Typical plasma density signals during propagation and 
after wave attenuation are shown in Figs. III-8a and III-8b. 


Diamagnetic probe measurements were made similar to those which have 

2 5 6 

been previously reported in studies on ion cyclotron wave heating. ' ' Because 
of the 360 Hz ripple inherent in the dc magnetic field power supply, the diamag- 
netic signal is passed through a high-pass filter before processing the oscillo- 
scope plug-in amplifier. It is then amplified and integrated and the resulting 
signal displayed, a typical display being given in Fig. Ill— 9 . 


A plot of diamagnetic signal versus B q at 0.5 mTorr is shown in Fig. 

Ill- 10. It is seen that r\%T L decreases gradually above 700G and then sharply 
down to 50G. The decrease is regular, and there appears to be no peak in 
near the wave resonances, nor dip where the stop bands occur. Figure III- 11 
depicts at 1 . 0 mTorr and corresponding values for electron density and 

magnetic probe signal. Although there is a close correspondence between density 
and b^, nX-"T^ gradually falls toward zero with no particular relation to the 
other curves. Idential behavior was registered throughout the entire 0.3 - 1.2 
mTorr range, and the only close correspondence was found at 0.3 mTorr, where 



Fig. Ill- 7 . Electron Density and b (p-p) at 0.8 mT. 
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a . 


P = 0.7 raT Z = 40 cm 

H 2 gas Oscilloscope Sweep: 200^isec/cm 

B = 400 G 
o 


b. 


p = 0.7 mT Z = 40 cm 

H 2 gas Oscilloscope Sweep: 200 ^isec/cm 

B = 50 G 
o 

Fig. Ill- 8 . Electron Density Signal from Microwave 
Interferometer during Propagation (a.) 
and Attenuation (b.) of Low Field Waves. 





p = 0 . 8 mT 


Z = 50 cm 


H gas Oscilloscope Sweep: 20 yusec/cm 

La 

B = 900 G 
o 


Fig. Ill- 9 . Photograph of Typical Diamagnetic 


Signal Trace. 



Magnetic Flux Density at Stix Coil (kG) 
Fig. Ill- 10. Diamagnetic Signal at 0.5 mT. 



39 


the wave phenomena was not exhibited. In this case, density, magnetic signal, 

and diamagnetic signal all decrease with decreasing magnetic field. 

• • 

It is interesting to examine the data by plotting curves of b^ versus 
neutral gas pressure for constant magnetic field B . Such a curve is shown in 
Fig. Ill- 12 for B q = 1.0 kG . The magnitude of magnetic signal gradually increases 
and peaks at about 0.8 mTorr, then gradually decreases. The bars in Fig. Ill- 12 
show the fluctuation of signal above 0.7 mTorr, which implies that although the 
amplitude of the propagating wave could be quite large , the wave tended to 
become cutoff more easily as pressure was increased. This behavior illustrates 
the fact, previously alluded to, that above about 1.0 mTorr there is a gradual 
"smearing out" of wave effects until all such effects are lost around 8-10 mTorr. 

The peak in b occurs in all the curves of b versus pressure. As B is 
z z o 

made lower, this peak gradually shifts slightly lower in pressure until it occurs 
at about 0.6 mTorr at 100G. At 300G, a double peak appears and a very strong 
cut-off region of neutral pressure occurs centered at 0.8 mTorr. This appears due 
to the widening of the cut-off regions discussed in Part A. The cut-off region 
is then displaced by the widening central propagation band and falls lower in 
the magnetic field range as the lower propagation band gradually disappears 
so that above 0.8 mTorr wave propagation once again occurs at 300G. This same 
effect is not as evident at other values of magnetic field between 0. 1 and 1 .0 kG. 

The peaking of the electron density versus pressure curve is also quite 

pronounced, as seen in Fig. Ill- 13, at B q = 100G. It is notable that this peak 

occurs about 0.2 mTorr higher than the corresponding peak in b^. In general, 

however, the agreement between the two curves is quite close, as seen in 

Fig. Ill— 14, at B = 500G, although there seems to be very little shift in the 
o * 

location of the peak density along the pressure axis. At 300G, the dip in b^ at 
0,8 mTorr mentioned above parallels a corresponding minimum in electron density, 
which also exhibits a double peak behavior. 
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C . Low Magnetic Field Harmonic Frequency Phenomena 

4 

Work on nonlinear plasma phenomena is extended here to the low magnetic 

field case. The method of detection of frequency spectra is identical to that 

4 

reported previously. The magnetic signal is terminated in a 50 ohm coaxial 
load, then passed through a 100:1 attenuator into a plug-in-amplifier unit in a 
Tektronix Model 555 oscilloscope. Previously, a Pentrix spectrum analyzer had 
been used as the plug-in unit, but a new Tektronix Model 1L20 spectrum analyzer 
has been obtained and was used in this series of experiments. 

A typical photograph is shown in Fig. Ill- 15a. The lower trace shows 
the magnetic signal being analyzed, while the upper trace gives the spectral 
content of the waveform. The far right vertical line of the graticle represents 
the lower end of the frequency scale, about 5 MHz, and the dispersion is 
5 MHz/cm. The upper limit of the detector, in this case, is about 52 MHz, 
about half a centimeter from the left end of the scale. The spectral lines are 
broadened considerably because of the very rapid sweep required; the main rf pulse 
is only about 1 msec long. Although this is a compromise in that a slower sweep 
is desirable, the data is satisfactory since there is very little overlapping of 
the frequences which occur. Also, the reduction in amplitude due to the rapid 
sweep could be amply compensated for by the gain control on the analyzer. 

A number of frequencies seem to be present in the magnetic signal over 
the entire range of the experiments. However, during wave attenuation, particu- 
larly at low fields, these harmonics become very substantial in amplitude. An 
example is shown in Fig. Ill- 15b for the condition of 0 . 6 mTorr pressure at 
B — 50G. The small spectral line on the far right represents the 5.8 MHz fun- 
damental of the main rf transmitter; the second and subsequent "spikes" 

proceeding to the left are at harmonics of the 5.8 MHz fundamental, second 

4 

through eighth. As noted in the last report, a Tektronix Model 191 constant 
amplitude since wave generator was used in the identification of the various 
spectral lines . 



a. Spectrum (Upper Trace) 

p = 0 . 7 mT 

H 2 gas 

B = 1000 G 
o 

Z = 50 cm 



b. Spectrum (Upper Trace) 

p = 0 . 6 mT 

H 2 gas 

B = 50 G 
o 

Z = 50 cm 


III- 15. (a.) Typical Spectrum from Plasma, 

(b.) Spectrum During Wave Attenuation. 
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The occurrence of the harmonics lines seemed to be sensitive to magnetic 

field in a similar fashion to the magnitude of b and electron density. Fig. Ill— 1 6 

z 

compares values of magnitudes of the 5.8 MHz fundamental and first harmonic 

(11.6 MHz) at 0.8 mTorr as a function of B . The curves seem nearly the inverse 

o 

of one another in form, the maximum of the first corresponding to the minimum of 
the second. 

The ratio of the magnitude of the 11.6 MHz signal to the 5.8 MHz signal 

may be taken as a measure of the amount of distortion present in the wave. The 

graph of this ratio at 0.8 mTorr is shown in Fig. Ill- 17. It is interesting to note 

that this curve is an inverse to the plots of b and density at 0.8 mTorr; i.e. , 

2 • 

the minimum distortion seems to occur in the range of 400-600G where both b 

z 

and density are at a maximum. Other harmonics of 5.8 MHz are also present, of 
course, but their behavior generally corresponds to that of the first harmonic. 

It therefore appears that maximum purity of spectrum occurs when wave propaga- 
tion is strongest and electron density is high. 

As a verification that the various frequency effects which have been 
noted are actually due to plasma-wave interactions, a current sensor was con- 
nected in the tuning circuit of the Stix coil. This sensor consisted of a small 
Rogowsky belt of about 100 turns wounds around a plexiglas form, connected to 
a BNC jack, and placed around the input strap to the Stix coil. 

The current envlope and corresponding spectrum are shown at 0.8 mTorr, 

100G in Fig. Ill- 18a. In Fig. Ill- 18b, b and its spectrum are shown (the narrow 

2 

pulse is the 40 kW preionizer and was not in the time sequence being analyzed). 
The 5.8 MHz "spike" of the coil current spectrum is extremely large (it is hard to 
see in the figure, since the line is deflected off scale) with only very small har- 
monic lines, the double-peaked pulse on the far left being a transient in the 
spectrum analyzer and not a true spectral line. In contrast, there is a great deal 

of harmonic content in the b waveform. It is apparent that the harmonic phe- 

z 

nomena occur in the plasma rather than in the transmitter or matching circuitry . 






Magnetic Flux Density at Stix Coil (kG) 

Fig. Ill- 17. Ratio of Amplitude of 11.6 MHz Signal to Amplitude 
of 5.8 MHz Signal at 0.8 mT, H gas. 




a. Coil Current Spectrum (Upper Trace) 


p = 0 . 8 mT 

H 2 gas 

B = 100 G 
o 
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b. b Spectrum (Upper Trace) 
z 

p = 0 . 8 mT 

H 2 gas 

B = 100 G 
o 

Z = 50 cm 


Fig. Ill- 18 . Comparison of Spectrum of b^ and 

Stix Coil Current During Wave Attenuation. 
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D. Discussion and Interpretation of Experimental Results 

As mentioned in the introduction to Section III, the wave phenomena 
previously observed below 1.0 kG have been referred to as the so-called "fast" 
or right-hand circularly polarized hydromagnetic wave. More careful investiga- 
tion, as described in Parts A-C of this section, reveals a complex behavior of 
the plasma which may allow as many as three bands of propagation in the range 
of 0.1 - 1.0 kG at a number of neutral pressures from 0.3 to 1.2 mTorr in hydro- 
gen gas . 


By examining the cold plasma dispersion relation for a cylindrical hydro- 

7 

gen plasma bounded by a' vacuum, one finds 


10 > ■574-5 ( I o? 

r 0 



(III- 3) 


for propagation of the "fast" hydromagnetic wave, where CO is the frequency of 
oscillation, jQj is the ion cyclotron frequency, r is the radius of the plasma, 
and n g the electron density. At to = (24T) 5.8 (10)^ (driving frequency in 
radians of the main rf transmitter), and taking the plasma radius to be 3.0 cm, 
we may write equation (III— 3) as, 

n e > 2.28 do) 7 B 0 2 m . 4) 

where B is the static magnetic confining field, expressed in Gauss. For 

° 12 
B q = 500G, n^ must be greater than 5.7 x 10 particles per cubic centimeter, 

a density higher than any observed in the data taken at 500G. Thus, it would 

appear that only the very lowest propagation band (50 to 75 G) observed could be 

the "fast" wave. 

It should be emphasized that the laboratory plasma under consideration 

does not always fall under the classification of "cold" plasma. Moreover, radial 

density gradients may be present such that the peak density in the central core of 

12 -3 

the plasma is considerably above 5.7 x 10 cm , in which case "fast" wave 
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propagation would be possible along the central high density column. Thus, the 
propagation of the "fast" wave is not ruled out, although it seems unlikely at 
the higher values of B q (i.e. , for 300G). 

2 3 

A number of sources ' have discussed plasma heating by the "fast" wave. 

2 

Rothman and Chung report better than 70 eV increase in the (T + ) 

sum, although the theory is not clear on heating with the right-handed wave. In 

particular, the above reference includes diamagnetic probe data supporting the 

findings quoted. The diamagnetic probe and supporting equipment used in the 

5 7 

experiment at this facility are identical to those used previously ' to report 
plasma heating by the ion cyclotron wave. The data obtained was satisfactory, 
and the behavior of the signal was as expected near B = B c> Since the identical 
equipment was employed to take the data reported in Part B of Sec. Ill, the data 
appears to be authentic and to indicate a cooling of the plasma at low magnetic 
fields, which would not be characteristic of the fast wave. 

Since the findings in the three preceding paragraphs are in agreement, it 
is tentatively concluded that the wave phenomena observed are not predominantly 
the "fast" hydromagnetic wave. It is anticipated that use of the low-current 
dc power supply discussed in Part E will facilitate the verification of the above 
conclusion by providing accurate values of low magnetic fields. 

A number of possibilities have still to be considered among waves which 

might propagate at low magnetic fields . Among these are electron and ion 

acoustic waves and electrostatic ion and electron cyclotron waves. It should be 

pointed out that although these waves are primarily electrostatic in nature, there 

are small magnetic fields associated with them which could be detected using 

magnetic probes such as the b probe used here. 

z 

The most probable choices to consider are the ion acoustic wave and 

electrostatic ion cyclotron wave. Although no study of these waves has been 

attempted here for the case of a cylindrical plasma, the cold plasma modes for 

8 

the unbounded case are well known in the literature. Examination of the disper- 
sion relation for the electrostatic ion cyclotron wave reveals that for electron 
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temperatures of interest (e.g. ,^10 eV) , the mode propagates only for 3.8 kG, 

high above the magnetic field region of interest. 

The ion acoustic wave is found to propagate at some angle 9 to the 
magnetic field lines, or parallel to B q/ but perpendicular to B q only near the limit 
as B q — ► 0. However, at the driving frequency of the main transmitter, the 
wavelengths involved are considerably larger than the dimensions of the laboratory 
plasma so that any ion acoustic oscillation would be strongly influenced by the 
boundary conditions imposed on the plasma. 

Until more data is taken and theoretical work completed, the exact iden- 
tity of the low- field wave phenomena remains unknown. What has become clear 
is that such phenomena seems a great deal more complex and of considerable 
more variety than was previously thought. 

Reference to Fig. Ill- 14 in Part B reveals an optimum pressure such that 
electron density is maximized around 0.8 mTorr, at B q = 500 G. Since the elec- 
tron density gradually decreases above 0 . 8 mTorr (and this decrease was noted 
at all magnetic fields in the 0.9 - 1.2 mTorr range), the percent ionization 
decreases rapidly. Thus, it is seen that the total number of neutral particles 
will be three orders of magnitude or more greater than the electron density by 
10 mTorr. It is a reasonable result that in this pressure region, where ion- neutral 
collisions would be expected to become important, the low- field wave phenomena 
gradually disappear. It might be noted that around 8-10 mTorr the same "washing 
out" of ion cyclotron wave effects also occurs. 

Part C deals with the nonlinear frequency phenomena observed at low 
magnetic fields. The amplitudes of most of the harmonic frequencies seems to 
increase markedly during attenuation of wave effects and diminish during strong 
propagation. This would seem to indicate that attenuation of the propagating 
wave is due to damping, either collisional or collisionless, from plasma-wave 
interaction (i.e. , that plasma and wave are interacting). Because of diamagnetic 
signal measurements, however, from which a net temperature decrease is indicated, 
it appears that damping phenomena are not responsible for the wave attenuation. 
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Apparently wave cut-off occurs when propagation is no longer allowed by the 
varying parameters (n» B q ) and the boundary conditions imposed on the plasma, 
and the wave evanesces. The harmonic frequencies which then occur must be 
attributed to interaction of evanescing fields with the decaying plasma. 

E. Plans for the Next Report Period 

In the lower portion of its output range the 300 kW dc power supply is 
difficult to adjust to within + 50 Gauss . The present plans call for substituting 
a low-current (0-100 a) dc power supply which will facilitate more accurate and 
easily controllable measurements in the region of 0-1000 G. Some of the measure- 
ments discussed in Parts A, B, and C will be repeated, and the series of 
experiments will be extended to deuterium. The magnetic field range in the 
experiments will be raispd to include the top of the third apparent propagation 
band. 

It is anticipated that further modifications in the microwave equipment 
will eliminate the noise problem and the need for modulation of the klystron. 
Electron density, diamagnetic signal, and harmonic frequency data will be taken 
as described above. 

Investigation of applicable theory to the low- field propagation will be 
undertaken, and study of pertinent nonlinear theory leading to a better under- 
standing of the harmonic and frequency phenomena will be continued. Spectrum 
data will also be taken in the region of propagation and damping of the fundamen- 
tal ion cyclotron wave . 
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